Electromagnetic and phonon modes for superfluid He^ with a disk resonator 



O 
(N 



oo 



5-H 

CD 



a 



X3 
O 

o 



> 

OO 
(N 
in 

vn 

o 
o 



X 



V.M. Lokte-Jil and M.D. Tomchenkclil 
Bogolyubov Institute for Theoretical Physics, 14-b Metrologicheskaya Street, Kiev, 03680, Ukraine 

(Dated: May 31, 2010) 

We find the distribution of the electromagnetic field inside and outside a dielectric disk resonator 
placed in He-11. It is shown that this field consists of a collection of "circular" (c-) photons. The wave 
function '^c of a c-phonon for the He-II + disk system is calculated in the zero-order approximation 
in interaction. Due to the symmetry of the problem, the structure of ^c is such that a c-phonon 
possesses, similarly to a c-photon of the resonator, a definite energy and an angular momentum with 
respect to the disk axis, but it does not possess a definite momentum in the disk plane. 

PACS numbers: 07.57.-c; 71.10.-w 



I. INRODUCTION 

In the recent years, some interesting and, in a cer- 
tain sense, unexpected effects were discovered in the 
experimentsii^. Namely, a supernarrow absorption line 
at the frequency of the roton minimum was registered in 
the spectrum of a dielectric disk SHF resonator placed 
in liquid Hc-4. In an external constant electric field, the 
line is split into two ones. At the switching-on of a heat 
gun directed along a tangent to the disk, the absorption 
line is transformed in an emission line. These effects have 
no explanation yet, though the line itself is related, un- 
doubtedly, to a high density of states of He-II at energies 
close to the roton onc^. 

To explain the origin of the line and its specific fea- 
tures, it is necessary to determine, first of all, the elec- 
tromagnetic (EM) field of a resonator, as well as the wave 
functions of a phonon and a roton for helium with an im- 
mersed disk. The present work is devoted to this prob- 
lem. 



II. ELECTROMAGNETIC FIELD OF A DISK 
RESONATOR 

In the experiments described iniiS, a variable inhomo- 
geneous field with the strength i? < 10^ V/m was induced 
in a resonator. The field was mainly concentrated in a 
disk and created the deformations of a resonator which 
are pulsating in time and space. However, at the attained 
values of E, the total deformation of a disk was very small 
— at most O.lAi^ for a quartz resonator. Similar weak 
pulsations should play no role in the phenomena under 
study. Therefore, it is obvious that a roton is excited 
by the SHF field of a circular EM wave pulsating on the 
rim of a disk, rather than by deformations of the disk. In 
what follow, we will calculate the EM field of a resonator. 

In the experiments, the sizes of disk resonators were 
approximately identical. In^ and^, the resonators were 
fabricated of quartz and leucosapphire, respectively. The 
results obtained for the shape and the width of a roton 
line are close, but the numbers of the azimuth mode (for 
the roton frequency) are different. Below, we will obtain 
the general formulas for the EM field of a resonator and 



analyze the numerical values for the experimental con- 
ditions in^. Let us consider the EM wave propagating 
in a quartz resonator with the shape of a disk with the 
thickness hd = \ mm and the radius Rd = 9.5 mm. The 
dielectric permittivity tensor Ejk for the quartz under 
study is diagonal in the coordinate system (CS), whose 
Z axis coincides with the geometric axis (it is also opti- 
cal) of a resonator; in this case, e^ — 4.63, and e±_ — 4.43 
in perpendicular directions^. 

In calculations of the EM field, we are based on the 
Maxwell equations in a medium: 



divT) = 0, dz-wB = 0, 



rotE 



laB 
Vdt' 



rotH 



19D 



D = £E, 



(1) 



(2) 



(3) 



where c is the light velocity in vacuum. For quartz and 
helium, /x ss 1, therefore, B = H. 

We now find the vector potential A connected with E 
and H by the relations 



E: 



19A 
B = TotK. 



Vipe 



(4) 



(5) 



We use the transverse calibration ipei = and pass into 
a cylindrical CS (CCS) p, ip, z with the origin at the 
disk center and the Z axis coinciding with the axis of 
a resonator. In the CCS, the tensor e^fe is diagonal: 
Exx — £yy — £±, £zz — Ez- For the field in quartz, rela- 
tions H]) and (gl) yield 



dwA(r,i) = ( 1 

^ _L 



dz 



fd{r), 



(6) 



where fd is some function independent of t. Since we are 
interested in EM waves, we can take fd = 0. With the 
help of ^ and ([5]), we obtain the following equation for 
A: 



AA- 



C2 9t2 



^-1 



, dA, 
dz 



= 0. 



(7) 



For quartz, the values of e^ and £\_ are close. Therefore, 
we can neglect their difference and consider that e^ — 
£± — Ed ~ 4.63, which simplifies the equation: 



AA- 



1 d^A 



0, ci_ 



\/eI 



(8) 



This equation has a solution A directed identically at 
all points of a resonator and another solution directed 
according to the symmetry of the disk with the p-, ip-, 
and z-components. It is natural to expect that a res- 
onator amplifies maximally those components of the field 
which correspond to its symmetry. It follows from the 
experiment^ that this is true, and, in addition, the prin- 
cipal components of the field E near a resonator are the p- 
and (^-components, whereas the value of the z-component 
is less by three orders. Therefore, we neglect the latter 
and consider that the field A in a resonator and in helium 
has only p- and t/j-components. 

The equation for the field outside a resonator (in he- 
lium) has the form 



AA- 



1 d^A 



dt^ 



= 0, 



(9) 



where Ch = c/ ^Jsh, Eh — 1.057 (here and below, h and 
d mean, respectively, helium and a disk). In order to 
determine A, it is necessary to solve Eqs. ^ and © 
with regard for boundary conditions (BCs) on the surface 
of a resonator: 



^\\ = ^11 ' 



-^11 "= K\ 



and, if there are no extrinsic charges. 



Bi = B'l, 



D"^ =D^, 



(10) 



(11) 



(here, the symbols || and _L indicate the relations to the 
surface, whereas the symbol _L in the other cases means 
the relation to the disk axis (the Z axis)). 

We now calculate the field A outside and inside a disk. 
The general form of a solution A(p, Lp, z) is unknown else 
and, generally speaking, complicated. In principle, the 
field can depend on the shapes and the sizes of a container 
and antennas'* (for example, in the experiments ini—, two 
antennas are positioned in the disk plane on two sides 
from it at a distance of 13 mm, i.e., ss 1.37i?d from the 
disk axis). The geometry of a resonator is such that the 
field inside a disk can be determined with the use of the 
separation of variables: 



A{p,ip,z,t) 



dQzS-QAPi^) 



X [bciQz] cos {Q,z) + bs{Qz) sin (Q^z)]] + c.c. (12) 

Here, we took into account that the observed field is 
real. The solution contains no sines, because the sys- 
tem is symmetric relative to the reflection 2; — > — z. It 



is known from the experiment that a disk enhances the 
field E mainly inside itself. Outside the disk, the field is 
slight, rapidly decreases, and disappears practically at a 
distance of 2 mm from the disk. Therefore, we assume 
that the structure of the solution outside the disk is such 
that we can approximately separate the variables z and 
(on the other hand) p, (p according to (H^). 

Since the z-component A is small, we can write 
aQ^{p,tp) = aQ^{p,ip)ep + aQ^{p,ip)e^. Then relation 
([8]) yields the equation for agj 



Ap,^aQ,. + iQ^YsLQ^ = 



(13) 



where Q'^ depends on Q^: 



Qi^^JiQ^Y-Ql g' = ^/ci = 2Wc±. (14) 

After simple calculations, we get the general solution of 
Eq. dUl) for a real Q'^: 



^qAp^v) ^J2j'^Q-9i{Qz)e'"' 



(15) 



X {[Ji^iiQ^p) - Q(g,)Ji+i(o»]ep 

+ ze^[Jz_i(Q» -I- ci{Qz)Ji+i{Qtp)]} ■ 

Here, I is an integer, Ji{x) is the Bessel function, and 
9i{Qz) and ci{Qz) are constants. The second indepen- 
dent solution of Eq. ([T^ proportional to the Neumann 
functions Ni±i{Qpp) is omitted, because it tends to in- 
finity as p — > 0. The radial wave number Q^ in (|T5]) is 
determined, according to ((HI), by the value of Qz', Qi is 
positive at Qz < Q'^ and imaginary at Qz > Q'^. For the 
imaginary argument, we have Ji{ix) = i^Ii{x^. The plot 
of the function J; (a;) for / = 67 is given in Fig. 1. The 
function Ji{x) oscillates outside the disk, whereas Ii{x) 
increases monotonously and rapidly for all x. The EM 
pumping field creates some field A with a given frequency 
v in the disk and outside it, and this field increases in a 
resonance manner at definite values of v, Qz, and I. We 
do not calculate the exact condition for a resonance and 
the width of the resonant Z-mode, because it is easy to es- 
tablish which modes of ([T5|) are observed with the help of 
experimental data. It can be expected that the approxi- 
mate condition for a resonance consists in the proximity 
of the EM field on the surface of a resonator to zero (see 

m)- 

Experiments revealed various resonance modes, from 
which the first z- and the first radial harmonics were 
studied in detail. By P; {rip — 1, 2, 3, . . .), we denote 
values of x, for which Ji {x) — 0. The first radial harmonic 
[up = 1) means that the field in the disk is distributed 
over the radius so that it includes only the first half- wave 
of the function Ji{x) and becomes almost zero near the 

disk edge. Therefore, Q'iRd ~ Pi ■ To be more exact, as 
p increases from zero to Rd, the field |A| increases firstly, 
attains a maximum near the disk edge (p ~ Rd — Q.?) mm). 
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Fig. 1: Solid line — the Bessel function J;,,^j+i(a;) = 

^67(2:); 000 — the Neumann function —Nq7{x) which 
grows very rapidly at x < 60. Therefore, the function 
1 9 {1 9 {~^ 67 {x))) denoted by • • • is shown in this region, 
the functions a/ J|g {x) + 



+ + + 



^qJIt{x) at s < a;o and 
at a; > a;o present 



the radial behavior of the field A for z = (in this case, 
xo = QiRdP^ 72.1 — disk edge, Qi/Qi ~ 0.34). 



then decreases, and takes ~ 1/8 of the maximum value at 
the disk edge {p = Rd)- The first z-harmonic means that 
the distribution of the field over z is close to cos [zir/hd)- 
On the upper and lower surfaces of the disk z = ±hd/2, 
the field A is less approximately by 8 times than that in 
the middle plane of the disk (2 = 0) at the same p and ip, 
which yields hd ~ l.0&7hd and Qz — t^ /hd- This allows 
us to uniquely determine Q^: Q^ = Qi w 72.1A'i/Rd- 
The BCs ^ and ^ yield Ap{p = Rd) = 0, which 
gives 



ci{Q,) = Ji-i{Q'iRd)/Ji+i{QtRd) 



(16) 



The roton line is observed for a certain mode char- 
acterized by Irot- In^, the quantity Irot was defined as 
the ratio of the frequency to the step between modes, 
and its value was estimated as Irot ~ 78. However, 
the approximate condition of resonance (J19p implies that 
the connection between the frequency v and I is not 
strictly linear, and the step must somewhat increase 
with decrease in v. By averaging |A| over time, we ob- 

tain |A| ^ fiQ.p) = ^J^^iQip) + cfJ^_^,iQip). Ac- 

cording to experiments, the maximum value /,„ of the 
function / on the interval p = ^ Rd is attained at 
Pm « i?d — 0.5 mm. For hd obtained above, such a value 
of pm follows from calculations for Irot ~ 66. In this case, 
fm « 0.17, and / sa 0.05 at the disk edge. The experi- 
ment gives that A at the edge is equal to 1/6-^1/10 of the 
maximum value at the same height, i.e. / sa 0.02 -^ 0.03. 
We now obtain the final solution for the field A inside 



the disk: 

Adip,(p,z,t) 



Ar 



^ilip—iujt 



cos (irz/hd) 



X {[Ji-iiQip) ~ ciJi+i{Qip)]ep (17) 

-I- ie^[Ji-iiQip) + ciJi+i{Qip)]} + c.c. 



Am — 



EmC 2MEmC 



2u}fm 

66, Q 



UJ 

0.265, 



Qi « 72.U3/Rd, 
hd - 1.087hd. (18) 



The experiment^ indicates that two counter circular 
EM waves propagate on the disk, and the amplitude of 
one of the waves is larger by 2 orders than that of the 
second one. Below, we will neglect the weaker wave char- 
acterized by a different sign of I. 

It is worth to note that the circular EM wave ([T7|) has 
no definite z-momentum. Indeed, acting by the operator 
Pz — —ihd/dz on state ([T7]), we obtain ~ s'm{TTz/hd), 
i.e. the state is changed. 

Thus, the system of waves in a resonator is charac- 
terized by three quantum numbers: Z, Up, and Qz- For 
the field in the disk, we possess solution PTI) . p^ with 
Qz — T^ /hd, Up = 1, and different I. Since the field is in- 
significant near the disk edge, the resonance frequencies 
v{l,np,Qz) are determined approximately by the equa- 
tion A{p = Rd) = 0. Whence we get Ji_i{Q'^Rd) = 0, 



(rip) 



i.e. Q^pRd ^ fi\Z' , or 



'2TriyRd 



£± - (QzRd) 



2 _ ,,("p) 
t^i-i 



(19) 



This is an approximate condition of resonance. For each 
mode v{l,np,Qz) (HH), the distribution of the field A in 
the disk at large I is similar to a circular gallery. Such 
resonance modes are called "whispering-gallery modes" , 
because it was noted else in ancient times that a word 
pronounced by whisper at certain places of a circular 
gallery at large temples is heard at a remote part of a 
temple. 

We are interesting in the modes v^l^l^n/hd)- At 
T=1.4K, the roton absorption line was observed at the 
frequency Vrot ~ 180.3 GHz which corresponds, as shown 
above, to I — Irot ~ 66. We consider that, at Z ^ 1, the 
following relation is true^: 



M/ 



\^^ w/ + 1.856^1/^ + 1. 033 r^/^ 



(20) 



(1) 



/^i^a^ 



72.719. Let us denote the 



This yields p\ 

harmonic v(lrot,^,''^ /hd) by vi^^^. For the sizes Rd and 
hd taken fromi, relation pj)) yields vi^^^ ~ 1.007i^,-ot = 
Vrot + 1.26 GHz. In the limits of the roton azimuthal 
mode, the frequencies differ from Vrot by ±1.2 MHz, 
whereas the frequencies vi^^^ for a resonator in He-II and 
in vacuum differ by ~ 20MHz^. Therefore, the exact 
condition of resonance must give vi^^_^ — Vrot i 1 MHz. It 
is easy to prove that condition (|19l) is sufficiently close to 
the exact one. 



Consider the field A in helium. Near the disk, it sat- 
isfies Eq. (|ni), whose solution at ^2 = looks as 



Ah{p,Lp,Z,t) 



E 

I 



dQMQz)F{Q,,z)e 



ilip—iujt 



X {[Ji_i(Q»+aKQ.)^i-i(Q»](ep + *e^) 

+ [£i{Q^)Ji+i{Q';p) + bi{Q,)Ni+,{Q';p)] (21) 

X (-Gp + ie^)} + C.C. 

To determine Ah, we use solution (fT7| . ([TSl) for the field 
inside the disk and the BCs dTUl) and ([TT|) . 

a) Regions above and under the disk, \z\ > hd/2,p < 
Rd- Here, we neglect the Neumann functions in ([Si]) , be- 
cause they increase unboundedly as p — >■ 0. As solutions 
for F{Qz,z), we can take functions of the form e*"^ or 
e"^. Relations (|10l) and ([TT|) imply that the solutions A 
on the disk surface must coincide for the disk and helium, 
A/i = Ad- Therefore, relation (PT|) is reduced to the form 

Aft = 22ILe^ii'-p-'^t)^-K,(\z\-hd/2) 
8 
X {[Ji^iiQ'lp)~ciJi+i{Q'lp)]ep (22) 

+ ie^[Ji-i{Q'l p) + ciJi+i{Q'lp)]} + c.c. 

The sewing A/i = A^ on the disk surface requires that 
Q^ coincide with Qi from IH]). For (gH, we have Qf = 



-E/i + K^, and the condition Q^ = Qi gives k^ ~ 

6.53//id. 

b) Region in helium around the disk, \z\ < hd/2,p > 
Rd. Relations ^ and ^ yield A^ = A^ and AJJ = 
A = 0. Then only the harmonic with Q^ — Tr/hd and 
/ = 66 remains in ()2ip in the sum ^ J dQ^ , and the 

function F(Qz, z) is reduced to cos (zTr/hd). In this case, 
we have for the roton frequency: 



Qi = ^JiQ'^y ~Ql~ 24.626/ Rd, 



Q'' = uj/ch, Qz = Tr/hd. 



(23) 



(24) 



At such Qi, values of the Neumann functions (see Fig. 
1) in (prt are greater by 20 orders than values of the 
Bessel functions for the region with helium near the disk. 
Therefore, the latter must be neglected (by the physical 
reasoning, solution ([21]) should be written in terms of the 
Hankel functions; since the Bessel functions are small in 
them, only the Neumann functions remain). As a result, 
we obtain 

Aft « g-AnC'''^-'^' cos (TTz/hd) 

X {[Ni^i{Q'lp)-biNi+i{Q'lp)]ep (25) 

+ ^e^[iV^_l(g'lV) + biNi+iiQ^^p)]} + c.c. 

The functions Ni{x) have the asymptotical 

N,ix^O)^-^J—^(-) ^Nr{x). (26) 

TT \X J 



The numerical analysis indicates that, for values of 
the argument x ~ QiRd ~ 24.6, this asymptotics is 
approximately (with a correction coefficient) satisfied, 
namely: Nfi^{Q>lRd) « 10.38iVg"f (Q^'ii'd) « -1.59 • lO^o, 
NeriQ'lRd + 1mm) « 17.65N^j'{Q'lRd + 1mm) « -3.3 • 
10", NesiQ'lRd) ~ U.2Nii{Q'lRd) « -6.05 • lO^s, 
A^65((3ii?d + lmm) w 19.39Ngi{Q'lRd + lnim) « -1.57- 
10^'^. The condition A^ == A^ = is satisfied at hi == 

Ni^i{Q'lRd)/N,+i{Q'lRd), which gives bge « 0.0381. We 
can avoid great numbers in solution (j25|) . if it is rewrit- 
ten in the normalized form and by taking the condition 
A = A^ into account: 

Aft « AnJi-i(Qii?d)e'"^~'"*cos(^z//i<i) 

X {K_i(Q^p)-6,nj+i(Q^p)]ep (27) 

+ ie^K-i(Qip) + 6;n/+i(Q5'p)]} + c.c, 

where I = 66, Ji-i{QiRd) ~ 1/27.831, and 

"/±i(QiP) = Ni±iiQ'lp)/Ni^iiQ'lRd). (28) 

For the region with helium, relation (I27p can be approx- 
imately written near the disk in a simple form 

Aft « ie^A„2J,_i(Qii?rf)e^''^-^"* 

X cos {ttz / hd){Rd/ pY^^ + c.c. (29) 

c) In the region \z\ > hd/2,p > R, we sew together 
the solutions for the regions \z\ > hd/2,p < R and 
1-^1 < hd/2,p > R along the surface of their intersec- 
tion. This surface is symmetric relative to a turn around 
the z axis and intersects any of the planes z, p along a 
certain curve z(p) which cannot be calculated analyti- 
cally. Moreover, the analysis indicates that the intersec- 
tion happens not for all z and p. This means that the 
solution is more complicated in this transient region and 
cannot be determined by the separation of variables. Be- 
low, we will use a rough sewing, by considering that there 
exists a line z(p), along which a smooth sewing is realized. 
Such an approximation is apparently admissible, because 
the field is small in this region. The final solution for the 
field in helium near the disk has the form 



Aft = A^ 



(30) 



A(J « ^„e'("^-"*)[ai(p,z)(ep + ze^) 

+ a2ip, z){-ep + ie^)] , (31) 



ai{p,z) 



a2{p,z) 



\e-^'^\^\-^-'^J Ji^^{Q,p) (/), 
^cos{zT:/hd)ni_i{Q\p) {II), 



27.83 



^^^cos{zTT/hd)ni+i{Q\p) {II), 



(32) 



(33) 



where /,// stand for the regions (/ : \z\ > hd/2,p = 
-=- p{z); II : p> Rd, |z| = -=- \z{p)\), and z{p) or p{z) 



is the sewing line. In this case, / = 66, cgg w 0.265, 
bee « 0.0381, ha « 1.087/irf, k^ « 6.53/hd, 1/27.831 = 
Ji-iiQiRd), Qi ~ 72.143/i?d, Q^' « 24.626/i?rf, A^ = 
2.94:Emc/u}, Em — 10* V/m (value of Em in the Interna- 
tional System of units is taken from the experiment* for 
the frequency band of a generator AVpump — 50 kHz). 

The presented solution is in an approximate agreement 
with experiment. Only one difference can be noticed: ac- 
cording to the experiment, the field Ah decreases by 1-2 
orders, as the distance from the disk increases by 1 mm. 
From ([50)1 - (|33 p . we obtain that the attenuation in regions 
I and II is as high as ~ 700 and ~ 43 times, respectively. 
That is, the attenuation is too strong in region I. How- 
ever, since we used the approximate experimental data 
on the field, the divergence can be related just to this 
circumstance. 

In practice, each resonance mode is a very narrow band 
of frequencies, for which the EM field has shape of a 
"dome" . Most likely, this testifies to the resonance am- 
plification of modes with some dispersion of Qz and Qp 
(independently). Respectively, the resonance frequency 
is eroded with the formation of a dome. But solution 
([50)) - (|55|) does not consider the dome and implies that 
the EM field has a single frequency, so that the coeffi- 
cients are found for the frequency I'rot — 180.3 GHz. 

According to quantum electrodynamics-, in order to 
quantize the electromagnetic field, it is necessary to know 
the photon wave functions (WF) "^p^ot ("- i^ ^ collection 
of quantum numbers characterizing a state of a photon) 
which form the basis, in which the general solution of the 
Maxwell equations for a specific physical system is ex- 
panded. The collection of basis functions depends on the 
symmetry of the system. Therefore, photons are of dif- 
ferent types — plane, circular (or cylindrical), or spher- 
ical — and are characterized by different collections of 
quantum numbers. If the system is translationally in- 
variant, then it is convenient to expand the field A in 
plane waves. In this case, a photon has a definite momen- 
tum ftQ and a definite energy fko, and Q — uj/c. In the 
case under consideration, the disk violates the transla- 
tional symmetry. However, the circular symmetry holds. 
Respectively, a solution of the Maxwell equations for A 
takes form ^7^, dUl), ([30]) -([331). Whence we determine 
the WF of a photon with Qz — 7T/hd,np — 1 for the p- 
and (^-polarizations: 

*pfot = e*'^"*"*cos(zVM[Jz_i(Qip)TQj/+i(Qip)] 

(34) 
inside the disk, and 



phot 



s''^-»-*(ai(p,z)Ta2(p,^)) 



(35) 



ouside the disk, where the upper sign in the parentheses 
is related to the p-polarization (normalizing factors are 
omitted). 

With regard for the angular momentum operator L = 
—ih[r X V] (in particular, Lz = —ihd/dip), it is easy to 
prove that the WF '^phot is characterized by eigenvalues 
E = huj and Lz — hi. However, the momentum for states 



(p4)l . ((35)1 is not defined. We call such states "circular 
photons" (c-photons). 

Hence, a resonator creates some number of identical 
c-photons with I = 66, Up — 1, and Qz = T^/hd- In 
space, a c-photon is localized in the disk and near it. 
We note that such a photon cannot be represented as 
a superposition of plane photons. Indeed, let the EM 
field A/i in helium be expanded in plane waves with the 
wave vector k — Chjio. Since the disk and helium have 
different values of e, a photon, being plane in helium, 
is reflected from the cylindrical surface of the disk in 
the form of a fan of diverging almost radial waves and 
is refracted in a complicated way inward the disk as a 
system of converging waves. Therefore, photons are not 
plane in helium or in the disk. As follows from the laws 
of conservation, a quasiparticle created by a c-photon 
in helium must have the same quantum numbers as the 
c-photon (energy and angular momentum), but it has 
no momentum in the disk plane. This implies that a 
phonon created by a c-photon in helium must also possess 
the circular symmetry. In this case, a c-photon emitted 
by a resonator can be approximately represented as a 
collection of radially moving almost plane photons, the 
last being wave packets with size ~ \phot- Such photons 
can create plane rotons or phonons, if the disk or, as was 
assumed in^, helium absorbs a recoil momentum. But 
this is already a combined process, and its probability is 
much less than that of the direct c-photon — > c-phonon 
process. 



III. CIRCULAR PHONONS 

In helium far from the disk, ordinary "plane" phonons 
and rotons, being wave packets localized in space, are 
propagating. But, near the disk and also far from it in 
the case where A of a phonon is of the order of the disk 
size, the structure of a phonon must correspond to the 
symmetry of the disk. 

As an example, we consider a free particle in the space 
with an infinite cylinder with radius Eic- We assume that 
the particle cannot penetrate into the cylinder. There- 
fore, its WF *(r,i) satisfies the BC 



*(p = i?e,t) = 
and the Schrodinger equation 

dt 2m 



(36) 



(37) 



In the stationary case where ^{v) = e*"*^(r,t), the 
Schrodinger equation takes the form of a wave equation 

A* + k'^^ = 0, fc2 = 2muj/h. (38) 

Solutions of Eqs. (p6l) and (p8|) are the functions 



^(/, kz, kp) = e-'^''e'^[aiH[^\kpp) + biH\''\kpp)], 



(39) 



where kf, = k"^ - kl, h\^\x) = J; (a;) + iNi{x) and 

HI {x) = Ji{x) — iNi{x) are the Hankel functions, and ai 

and bi are selected so that aiH^ '{kpRc) + biHi '{kpRc) — 

0. In view of the asymptotics Hi {x -^ oo) = \f2pKx ■ 

exp{ix — ml/2 — in /A) and H^ [x — > oo) — ^J^jwx ■ 
exp{—ix + i'Kl/2 + in/A), these functions describe the di- 
verging and converging waves, respectively. Thus, if an 
impenetrable cylinder is present in space, the solution 
for a free particle is represented by circular waves p9|) 
(with various I and k^), rather than plane ones. The so- 
lution differs from a plane wave, because the interaction 
is indirectly present through BC. 

If a disk is present instead of a cylinder, and ^I^ = on 
its whole surface, then the solutions of Eq. psp take the 
other form; 



A-2 



^{l,k„kp) 



/c , kz 



^ — ikz 



+ e''^'')e''^Ji{kpp), 



(40) 



27r 



0,±1,±2,.. 



(41) 

and Ji{kpRd) = 0. Like that in Section 2, the Neumann 
function Ni is not present in the solution, since it in- 
creases unboundedly near the disk axis. 

We now consider helium surrounding the disk. The 
microscopic model for He-II without disk is constructed 
in the main (see, e.g., survey^) for periodic BCs, as the 
volume of the system tends to infinity. The model in- 
volves the WFs of the ground state ^o and a state with 
one phonon ^c^o- A specific feature of our problem con- 
sists in the presence of a disk in helium. It would be most 
proper to find the functions ^o and ^c'J'o with zero BCs 
realized in the nature and with regard for a disk. This 
requires to construct the full microscopic model of He-II 
with a disk, which is a very complicated problem. There- 
fore, we limit ourselves by the calculation of ^c for an 
infinite system without regard for BCs. But, in this case, 
it will be necessary in certain situations to pass from J dk 
to the sum ^ and to know the value of ^c^o on bound- 

k 

aries. To his end, we will consider that, according to the 
preliminary analysis, the zero BCs lead to the equations 



sin(A:^|z| + a)\^^±h^/2,±H/2 =0, 



(42) 



MkpRd) = MkpR^) ^ 0, (43) 

which yield the conditions of quantization for kp and kz : 



kz = 



H-K, 



n, :=±1,±2 



, -i-i-, ..., 



kp — 



^oo — Rd 



>Z 



(44) 



(45) 



(n^ ^ 0, since ^c^o will not be zero on the z-boundaries 
otherwise). Here, Roc is the radius of a container with 



helium, and Rd depends on Up-. for the least Up — 1, re- 
lation PS|) yields Rd ~ l.5Rd, and Rd decreases to Rd 
with increase in Up. At small Up ( < I), there exists no 
solution fcp, for which the relations Ji{kpRd) = and 
Ji{kpRoo) = from (^5)) would be simultaneously satis- 
fied. However, the symmetry of the system deviates from 
the cylindrical one near the container walls. Therefore, 
the relation J;(fcpi?oo) = should not apparently hold, 
and only Ji{kpRd) = is valid. This yields 



kpih rip) = ^J.i'''/Rd, rip = 1, 2, 3, 



.. < I. 



(46) 



The last relation can be rewritten in the form of (|45p . 
by introducing Rd- Though conditions (|42|) -(j46 |) will be 
used, we will find the WF of a phonon in a simpler ap- 
proximation, by neglecting the BCs (as usually the mi- 
cromodels of He-II are constracted^.) . 

It follows from the A^-particle Schrodinger equation 
that if the WF of the ground state of He-II is represented 
in the form ^o — const ■ e^ , then the WF V^k of a plane 
(p-) or circular (c-) phonon satisfies the equation 



N 



2?77,4 



^[V^ + 2(V,5)V,]vI/k = E{k)^^. (47) 



We now consider the zero approximation without any 
interaction between atoms. In this case, S'^const, and 
(|Tfl) is reduced to 



N 



£-Ev,^^k = i.(k)*k 



2m4 



(48) 



which is the Schrodinger equation for N free particles. 
For a translationally invariant system in the case where 
a single particle has a momentum Kk, and the rest ones 
are immovable, a solution of the equation looks as 



*i 



P-k 



1 ^ 



,ikrj 



21.2 



Eo(k) 



2m4 ' 



(49) 



where p^ are the "plane" collective variables, and m^ is 
the helium atom mass. This solution serves as the zero 
approximation for the WF of a p-phonon. The consid- 
eration of the interaction leads, as known, to the trans- 
formation of the one-particle excitations P^ to collec- 
tive ones: ^k acquires corrections nonlinear in pk, and 
-Bo(k) ^ k^ is replaced by a more complicated dispersion 
law for quasiparticles. If a disk is present in helium, the 
translational symmetry is broken, but the circular sym- 
metry holds. Therefore, according to (pS]) and (|40l) . the 
solution of (SHI) is the WF 



*c(^,fc^,fcp 



CLk^,k 



N 



ip ^e*(i'^.-'^t)e'fe^^.J;(fcpp) 



j=i 



N 



Xl'^'.'=^-'^p('"j) = Pc{l,kz,kp). (50) 



It is the zero approximation to the WF of a circular 
phonon in helium-II with the immersed disk (the sum- 
mation is made over all atoms). We omit the Neumann 
function, since namely function (j50p is a solution under 
the most correct zero BCs. The consideration of the in- 
teraction between atoms will lead to the appearance of 
corrections to (j50p which are nonlinear in pc, but we will 
neglect them. In ([50]) . pc means the circular collective 
variables. 

The following question is of importance: Does the en- 
ergy -Ec(k) of a c-phonon coincide with the energy i5(k) 
of a p-phonon at the same kl For a free particle (Eq. 
dSH])), the energy depends only on k (but not on k^ and 
kp separately). Analogously, the energy of a c-phonon 
must depend only on k. But, at kp <^ k, a c-phonon is 
close to a p-phonon, and, hence, its energy must be close 
to the energy of a p-phonon, by differing proportionally 
to smallest kp/k. Therefore, for any other kp, the en- 
ergy of a c-phonon must also be close to the energy i?(k) 
of a p-phonon. Generally speaking, the exact equality 
Ec(k) = E{k) is possible as well. 

Acting by the momentum operator Lz = —ih'Y^ d/dipj 

j 
on the WF of a c-phonon (|50p . we verify that a c-phonon 

possesses the intrinsic momentum L^ — hi. Thus, in 
what follows, we will use the zero approximation (j50p 
and the conditions of quantization pi)) -(P5 )) for the WF 
of a c-phonon. 



shape. But the analytic calculation of a new shape is 
a quite difficult problem, and we will neglect the differ- 
ence of the symmetry of the system from the cylindrical 
one, by considering that the container with helium is a 
cylinder with radius Roc and height H (then the helium 
volume V = ttRI^H). 

As known, the pair correlation function (7(ri,r2) de- 
termining the probability for atom 1 to be at a point ri 
and for atom 2 to be at a point r2 is presented by the 
integral 



dvs, . . .dvN'^l = g{Yi,V2) /V"^ . 



(55) 



For a translationally invariant system, g(ri, r2) = g{Ti — 
r2). In our case, a disk positioned in He-II breaks the 
translational invariance. But, at small |ri — r2|, the func- 
tion ^(ri , r2) is determined by the interaction of the near- 
est atoms, so that it should depend in helium with a disk 
only on the difference ri — r2, if ri and r2 are not too 
close to the disk. For a region near the disk (at distances 
of about several interatomic ones), g{Yi,Y2) ^ g{Ti — Y2)- 
But it is a very thin layer which hardly influences the 
processes in bulk. Therefore, we accept that the relation 
g{ri,V2) — g{ri — r2) is always true, and, hence. 



dr.':! 



■ dTN'^l 



gi\r,-r2\)/V' 



(56) 



IV. NORMALIZATION OF THE WAVE 
FUNCTION OF A CIRCULAR PHONON 

As seen from ([50]) . we need to know the coefficient 
c;,fez,fep (below, c) for the WF of a c-phonon. We will 
determine it from the condition of normalization 



/ 



di7|*,*or = 1, 



(51) 



where dft = dri . . .dr^, and r^ are coordinates of the 
j-th atom. Using ([SD]) . we have 



dn\^,^or = h + i2, 



(52) 



/i = ^ j dmlY^jfikpp,) = c^ j dn^ljfikpp^), 

(53) 
h = ^Jdml2j2e^- 



JKVJl -^2 ) p^'=^ (^31 -^02 ) 



Jl<j2 



X Ji{kpPj^)Ji{kpPj^ 



(54) 



In the real experiment, the disk is positioned in helium 
between two long cylindrical rods-antennas located in the 
disk plane at a distance of 1.37 Rd from the disk cen- 
ter. Near the antennas, a phonon wave loses the circular 



For atoms which are not located at the disk surface, the 
relation 

S{k) = l + yJ[g{r)-l]e-^'^'^dr (57) 

obtained for translationally invariant systems is also 
valid. In addition, if ri is far from the disk, then the 
relation 

fdr2... dvN'^l = I dv2g{vi , rs) = l/V (58) 

is true. Indeed, this integral determines the probability 
to find atom 1 at the point ri, and it is obvious that the 
probability for all points far from the disk is the same. 
With regard for dSS]) . we obtain that integral (|55|) is 



h 



fioo 

y I dviJfikppi) = ^ pdpJfikpP) 







2c' 



TTkpJriQQ 



JDykprCoQ^ t J, 



(59) 



where 



B{x,l) = nx / ydyJi {y x). 



(60) 



We note that the Bessel functions satisfy the relation^ 

/ (mi) \ / (m2) \ 







R 



i?2 r I > n2 

-"i- ^1 f (rrai)N 



R 



J'M 



(61) 



where Smi,m2 is the Kronecker symbol, and Jl{x) — 
■j^Ji{x). First, we consider large kp. At kpp ^ Z, the 
condition (|45p is valid and the following asymptotic is 
true^: 



Ji{kpP) 



nkpp 



cos a\kpp) 



Af 



1 



8fcpp 



sin a [kpp) 



a\kpp) ^kpp-'-^-'-. 



■nl -K 

y 4 



(62) 
(63) 



In this case, the function Ji{kpp) performs many oscil- 
lations on the interval p = -^ Roo, and the following 
relation more general than ()6ip is approximately valid: 



pdpJi{kpp)Ji{qpp) « 5k^,q^^B(kpRoo,l)- (64) 



7rfc„ 



Moreover, at kpR^o ^ I, we have 
B{kpR^J)^^^^f^ 



Jl (KpRoc) 



(65) 



where i?oo is a value of p which is the closest to Roo and 
is such that kpRoo is equal to one of the zeros pf of 
the Bessel function Ji{x). According to (|43|) . ^oo = -Roo- 
Relation ([62]) yields 



Jlikpp) « -^ 



nkpp 



sm a 



(fc,p)] + O ((/cpp) 



-3/2 



For /9 = Roo, we have Ji{kpp) 
cosa^kpRoo) ~ 0, which yields 



(66) 
0; therefore. 



[JlikpR 



pJ-'-oo) 



\k„B.^:pl^ 



2/TrkpRoo, 



B{kpRoo,l) \k„R„ 



»i" 



II |fcpi?oo»/' 



25' 



;;2 



TTrCpTloo 



(67) 



(68) 



(69) 



At small kp {kpRoo < 0; (|46l) is satisfied. In this case, 
the function B{kpRoo, I) can be determined only numer- 
ically. The analysis indicates that kp from (|^5| satisfies 
the relation B{kpRoo,l = 66) > 0.9. In other words, for 
all kp, we may take B{kpRoo, I = 66) ~ 1. 



With regard for ((56)) and the well-known relations 






~ix cos ipiLiltp 



dip, 



Six) 



1 

2^ 



-, — ik: 



""dk. 



(70) 



(71) 



(/? is arbitrary; ([7^ follows from ([ST]) ). let us write the 
vector q in the CCS. Then integral ([54)) reads 



7ri?2 



dipqqpdqpPidpip2dp2[S{qp, -k^ 



X Ji{kppi)Ji(kpP2)Ji{qppi)Ji{qpP2)- 



(73) 



Using conditions ([45]) and ([46]) for qp, we pass from / dqp 
to the sum „ ^ i, J^- In this case, (?„ and fc„ in (1751) are 

yp 

quantized identically. With regard for (|64p , we finally get 

^ 252(5(fc)-l)B2(A:pi?„o,0 



nkp{Roo - Rd) 



(74) 



By using relations ([5Tj), ([52]), ([59]), and ([74]), we obtain 
the required result for the normalization of the WF of a 
c-phonon: 



C = C;,fe, fc„ 



O.bnkpRoo 



B + B^S{k)-l){l-Rd/Roo)-^ 



(75) 
where B = B{kpRoo,l)- This formula is approximately 
true also for arbitrary kp which is not quantized accord- 
ing to (H5]) and ([^5]) . We note that, while integrating, 
we do not consider that helium atoms cannot be present 
in the volume occupied by the disk, but taking this cir- 
cumstance into account does not practically change the 
integrals and result ([75]). 



V. CONCLUSION 

Thus, we have determined the distributions of the elec- 
tromagnetic field inside and outside a resonator, as well 
as the wave function of a circular phonon. Without these 
quantities, it is impossible to calculate the SHF absorp- 
tion spectrum of liquid helium which arises due to the 
creation of quasiparticles in helium by the field of a res- 
onator. In our opinion, just the mutual transformation 
of excitations with the circular symmetry (photon ^ 
phonon or photon <^ roton) allows one to understand 
the process of absorption in helium with an immersed 



disk resonator, when the momentum conservation law is 
formaUy broken, and it is necessary to determine which 
quantum numbers of created and disappeared quasiparti- 
cles must be conserved. The calculation of the probabili- 
ties of relevant transitions and the description of the phe- 
nomena discovered in experimental worka^*^ are planned 



to present in the subsequent publications. 
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